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Abstract. Notions of interpolating variational splines and Paley- Wiener spaces 
are introduced on a combinatorial graph G. Both of these definitions explore 
existence of a combinatorial Laplace operator on G. The existence and unique- 
ness of interpolating variational spline on a graph is shown. As an application 
of variational splines the paper presents a reconstruction algorithm of Paley- 
Wiener functions on graphs from their uniqueness sets. 



1. Introduction and Main Results 

The paper introduces variational splines and Paley- Wiener spaces on combi- 
natorial graphs. Variational splines are defined as minimizers of Sobolev norms 
which are introduced in terms of a combinatorial Laplace operator. It is shown 
that variational splines not only interpolate functions but also provide optimal ap- 
proximations to them. Paley- Wiener spaces on combinatorial graphs are defined 
by using spectral resolution of a combinatorial Laplace operator. The main result 
of the paper is a reconstruction algorithm of Paley- Wiener functions from their 
uniqueness sets using variational splines. 

The following is a summary of main notions and results. We consider finite 
or infinite and in this case countable connected graphs G = {V{G),E{G)), where 
V{G) is its set of vertices and E{G) is its set of edges. We consider only simple 
(no loops, no multiple edges) undirected unweighed graphs. A number of vertices 
adjacent to a vertex ti is called the degree of v and denoted by d{v). We assume 
that degrees of all vertices are bounded from above and we use notation 

d(G) — max d(v). 
vev(G) 

The space L2{G) is the Hilbert space of all complex- valued functions / : V{G) C 
with the following inner product 



veV{G) 



and the following norm 

11/11 = ll/llo 
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The discrete Laplace operator L is dGfined by the formula [4] 

where v ^ u nieans that w,m € V{G) are connected by an edge. It is known 
that the Laplace operator £ is a bounded operator in L2{G) which is self-adjoint 
and positive definite. Let (j{C) be the spectrum of a self-adjoint positive definite 
operator C in L2{G), then a{£) C [0, 2] . In what follows we will use the notations 

Wmin = inf CJ,CJmax= SUp UJ. 

For a fixed e > the Sobolev norm is introduced by the following formula 



(1-1) \\f\\t,e= {el + Cf/^f 



The Sobolev space Ht^e{G) is understood as the space of functions with the norm 
(jl.ip . Since the operator C is bounded all the spaces Ht^e{G), coincide as sets. 

Variational splines in spaces L2(R'^) are introduced as functions which minimize 
certain Sobolev norms [S^, [5]. Sobolev spaces in L2(M'^) can be defined as domains 
of powers of the Laplace operator A in L2{M.'^) |42]. To construct variational splines 
on a graph G we are going to use the same idea by replacing the classical Laplace 
operator A by the combinatorial Laplacian C in L2{G). 

For a given set of indices / (finite or infinite) the notation I2 will be used for 
the Hilbert space of all sequences of complex numbers y — G /, for which 

<00. 

Variational Problem 

Given a subset of vertices W = {w} C V{G), a sequence of complex numbers 
y — {Vw} G h,w e W, a positive i > 0, and an non-negative e > we consider the 
following variational problem: 

Find a function Y from the space L2{G) whicli has the following properties: 

1) Yiw)=y^,w&W, 

2) Y minimizes functional Y \\{el + Cy^^Y\\. 

Remark 1. It is convenient to have such a functional in the Variational Problem 
which is equivalent to a norm. Thus, if the operator C has a bounded inverse in 
L2{G) (it is a situation on homogeneous trees of order q + l,q > 2) then we will 
assume that e is zero. Otherwise we assume that e is a "small" positive number. 
In what follows we will write e > with understanding that e = 0, if the operator 
C is invertible in L2{G) and that e > 0, if £ is not invertible in L2{G). 

We show that the above variational problem has a unique solution Y^'^ . We 
say that Y^'^ is a variational spline of order t. It is also shown that every spline 
is a linear combination of fundamental solutions of the operator {el + £)* and in 
this sense it is a polyharmonic function with singularities. Namely it is shown that 
every spline satisfies the following equation 



at, 
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where {aw}wew ~ {ctw(Y^'^)}wew is a sequence from I2 and 6w is the Dirac 
measure at a vertex w e W. The set of ah such splines for a fixed W C V{G) and 
fixed t > 0,e>0, wih be denoted as y{W. t. e). 

A fundamental solution V{G), of the operator {el + Cy is the solution 

of the equation 



(1.2) (£/ + £)%» =5^, 



where 5yj is the Dirac measure at «; e V{G). 

It is shown in the paper that for every set of vertices W = {w}, every t > 0,e > 0, 
and for any given sequence y = {pw} € h, the solution Y^'^ of the Variational 
Problem has a representation 



where L^^^^ is the so called Lagrangian spline, i.e. it is a solution of the same 

Variational Problem with constraints lJ^"'(w) = dw,v,w € W, where 6w,v is the 
Kronecker delta. Another representation is 

where {aw{Y^'^)}w£W is a sequence in 12- 

Given a function / e L2{G) we will say that the spline Y^'^ interpolates / 
on W if Y^'^{w) = f{w) for all w G W. It is shown in the Theorem 2.4 that 
for a given function / € L2{G) its interpolating spline Y^'^ is always an optimal 
approximation (modulo given information). 

Remark 2. It is important to realize that for a fixed set W C V{G) and fixed 
t,e> 0, the correspondence 

(1-3) ^ {aUY,f")},y = {y^} G h, 

where Y^^'^ is a sphne, depends just on the geometry of G and W. In other words 
the map (1.3) is responsible for the connection between "analysis" on G and its 
geometry. 

Our main goal is to develop spline interpolation and approximation in the so- 
called Paley- Wiener spaces. 

Paley- Wiener spaces on M** are denoted PWt^(IR),w > 0, and contain functions 
/ e L2{M.) whose L2-Fourier transform 



has support in [— a;,w]. The classical sampling theorem says that if / € PWa;(M) 
then / is completely determined by its values at points k7r/u),k € Z, and can be 
reconstructed in a stable way from the samples /{kw/co) by using the so-called 
cardinal series 

s J, f kTi'\ sin(wa; - /ctt) 

/(^) = E/(- 
feez 



\ to J LOX — kw 
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where convergence is understood in the L2-sense. In papers [35], [36], [15], [18], 
[24] , [28] , sphnes were used as a tool for reconstruction of Paley- Wiener functions 
from their uniqueness sets. 

The Paley- Wiener spaces and in particular sampling problems in these spaces 
attracted attention of many mathematicians [I], [2], [6], [17], [22], [21], [19]. C. 
Shannon [38^ , IT , suggested to use the sampling theory in Paley- Wiener spaces as 
a theoretical foundation for practical problems in signal analysis and information 
theory. Since then the sampling theory found many other applications in particular 
in image reconstruction and learning theory [39j, [40] . Some of the ideas and 
methods of the sampling theory of Paley- Wiener functions were recently extended 
to the cases of Riemannian manifolds, groups, and quantum graphs [8], [9], [lOj . 
[TT] . [13], [H], [23]- [34]. Splines on manifolds and quantum graphs were developed 
in [32], [33]. 

To define Paley- Wiener spaces on combinatorial graphs we use the fact that the 
Laplace operator £ is a self-adjoint positive definite operator in the Hilbert space 
L2{G). According to the spectral theory ^ there exist a direct integral of Hilbert 
spaces X = J X{X)dm{X) and a unitary operator F from L2{G) onto X, which 
transforms domain of s > 0, onto Xg — {x E X\X'^x E X} with norm 

Mr)\\x^, - ( / A2^||a;(A)|||(,)dm(A) 

\Ja(C) 

and F{C f) = X'^{Ff ). We introduce the following notion of discrete Paley- Wiener 
spaces. 

Definition 1. Given an w > wc will say that a function / from L2{G) belongs 
to the Paley- Wiener space PWuj{G) if its "Fourier transform" Ff has support in 
[0,c.]. 

To be more consistent with the definition of the classical Paley- Wiener spaces we 
should consider the interval [0,w^] instead of [0,a;]. We prefer our choice because 
it makes formulas and notations simpler. 

Since the operator C is bounded every function from L2{G) belongs to a cer- 
tain Paley- Wiener space PWuj{G) for some oj G 'j{£) and we have the following 
stratification 

L2{G) = PM^<.._(G) = U PW^{G),PW^,{G) C PW^,{G),uji < uj2- 

Different properties of the spaces PWij(G) and in particular a generalization of 
the Paley- Wiener Theorem are collected in the Theorem 3.1. 

For a subset S C V{G) (finite or infinite) the notation L2{S) will denote the 
space of all functions from L2{G) with support in S: 

L2{S) ^ye L2{G), ^{v) ^Q,ve V{G)\S}. 

Definition 2. We say that a set of vertices U C V{G) is a uniqueness set for a 
space PWuj{G),u> > 0, if for any two functions from PWuj{G) the fact that they 
coincide on U implies that they coincide on V{G). 

Definition 3. We say that a set of vertices S C V{G) is a A-set if for any ip e L2{S) 
it admits a Poincare inequality with a constant A > 

11(^11 <A||/:^||,(^e L2(5), A >0. 
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The infimuni of all A > for which S* is a A-set will be called the Poincare constant 
of the set S and denoted by A{S). 

It is shown in the Theorem 3.4 that if a set S* C V{G) is a A-set, then its 
complement U — V{G)\S is a uniqueness set for any space PWu:{G) with uj < 1/A. 
Since L2{G) — PWt^^^^{G) every function in L2{G) belongs to a certain Paley- 
Wiener space and one cannot expect that non-trivial uniqueness sets there exist for 
functions from every Paley- Wiener subspace. But it is reasonable to expect that 
uniqueness sets exist for Paley- Wiener spaces PW^{G) with relatively small a; > 0. 
It will be shown (see Section 3) that for every graph G there exists a constant 
^Ig > 1 such that for < w < I^g functions from PWu:{G) can be determined by 
using their values only on certain subsets of vertices. Namely, it is shown that for 
any graph G spaces PWui{G) with 



have non-trivial uniqueness sets. A more detailed description of uniqueness sets 
will be given in a separate paper. 

The main result of the present article is obtained in Section 4 and can be stated 
in the following form. 

Theorem 1.1. 1) Assume that C is invertible in L^iG) . If S is a A-set then any 
f G PWui{G) with u) < 1/A can be reconstructed from its values on U — V{G) \ S 
as the following limit 

f = lim Yf'^fc = 2',ZeN, 

k^oo 

where Y^''^ is a spline interpolating f on the set U — V{G) \ S and the error 
estimate is 

<2j''\\f\\,-f = Alj <l,k^2\l gN. 

2) If the operator C in L2{G) is not invertible, then for any A-set S and any 
< e < 1/A, every function f G PWuj{G), where 

< CJ < - £, 

can be reconstructed from its values on U — V{G) \ S as the following limit 

f ^ lim Yl''/,k^ 2\lG N, 

where YjJ'J is a spline interpolating f on the set U = V{G) \ S and the error 
estimate is given by 

<27'=|l/|l,7 = A(c^ + £) < l,k = 2\lGN. 

We know two papers [12) . |16| . in which authors consider sampling on Z" and 
Zat and one paper [2^ were sampling on Z" was used to prove some deep results 
in discrete harmonic analysis. But our approach to the problem and our results are 
very different from the methods and results of these papers. 



UJ 

k,€ 
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2. Variational splines on combinatorial graphs 

We are going to use the same notations and the same Variational Problem which 
were defined in the Introduction. 

Theorem 2.1. For every set of vertices W = {w}, all k > 0,e > 0, and for any 

given sequence y = {Vw} G the Variational Problem has a unique solution. 

Proof. Consider the set Mo{W) C L2{G), of all functions from L2{G) whose re- 
striction to W = {w} is zero. This is a closed subspace of L2{W). 

Given a sequence of complex numbers y = {yw} G h, the linear manifold 
A4{W,y) of all functions / from L2{G) such that /(w) = is a shift of the 
closed subspace M.q{W), i.e. 

(2.1) M{W,y)=Mo{W)+g, 

where g is any function from L2{G) such that for all w G one has g{w) — y^. 

Consider the orthogonal projection ht^e of the function g E L2{G) from (|2.ip 
onto the space M.{){W) with respect to the inner product in Ht^e{G),t > 0: 

The function Y^^'^ ~ 9 ~ ^t.e is the solution to the above variational prob- 
lem. Indeed, it is clear that Y^'^ £ A4{W,y). To show that Y^'^ minimizes the 
functional 

Y ^ \\{el + Cy/^YW 

on the set A4{W, y) we note that any function from Ai{W, y) can be written in the 
form Y^'^ +ip, where V' G Mo{W). Since Y^'"^ = .9 — ^t,£ is orthogonal to Mf){W) 
in Ht.e{G) we obtain for any e C 

\\{eI+Cr'\YZ^y + an? = \\{eI+Cr'^Yr'yrMa\^m 

that means that the function Y^^'^ is the minimizer. 

The fact that the minimizer is unique follows from the well-known properties of 
Hilbert spaces. The proof is complete. □ 

The following result shows that every solution of the Variational Problem l)-2) 
should be a "polyharmonic function" with "singularities" on the set W . 

Theorem 2.2. For every set of vertices W — {w}, w € ^(G), every t > 0,£ > 0, 
and for any given sequence y = {yw\ G ^2, the solution Y^'^ of the Variational 
Problem satisfies the following equation 

(2.2) (e/-f£)*yj'^= ^ a„J^, 

where {aw\wew — {ctw{Y^e^y\w^w is a sequence from I2. Conversely, if a function 
satisfies equation \2. 2]) then it is a spline. 

Proof. If 5u} is a Dirac function concentrated at a point w £ W then for any 
(j) G L2{G) the function 

= - ^ 4){w)5nj 
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belongs to A4o{W) and because every solution of the above Variational Problem 
l)-2) is orthogonal to A^o(T^) in the Hilbert space H2t.,e{G) we obtain 



0= {eI + C)"^Y^,y{v){eI + cf^^{v). 

■uey(G) 

It implies that 

v£V(G) w^W 

In other words {el + Cf Y^'^ is a function of the form 



where a^iY^'^) — /Yj^'^,(5u,\ € l2- Thus we proved that every solution 

of the Variational Problem is a solution of (12.21). The converse is obvious. The 



Theorem is proved. □ 

A fundamental solution E21 G V{G), of the operator (el + £)*, is a solution 
of the equation 

(2.3) {eI + CYE^,,^,^S,, 

where 6^ is the Dirac measure at w G V{G)- The following Theorem explains the 
structure of splines and it follows from (|2.2p and linearity of the set of splines which 
is a consequence of the last Theorem. 

Theorem 2.3. For every set of vertices W = {w}, every t > 0, e > 0, and for any 

given sequence y = {yw\ G Z2, the solution Y^ ^'^ of the Variational Problem has a 
representation 

^t,e — 2^ Vw^t.,e 1 

wew 

where L^g^ is the so called Lagrangian spline, i.e. it is a solution of the same 

Variational Problem with constrains L^^^{v) — 5u],v,w G W, where Sw,v is the 
Kronecker delta. Another representation is 

(2.4) r.lT''' = E o.AY,'^s''')E^t,e, 

where {aw{Y^'y)}w€W is a sequence in I2. 

Now we are going to show that variational interpolating splines provide an opti- 
mal approximation. 

Definition 4. For the given 1^ C V{G),f E L2{G),t > 0,e > 0,/^ > 0, the 
notation Q{W, f, t, e, K) will be used for a set of all functions g in L2{G) such that 

1) 9{w) = f{w),w e W, 
and 

2) ||(e/ + £)*/2.g|| < K. 

It is easy to verify that every set Q{W, f, k, s, K) is convex, bounded, and closed. 
The next Theorem shows that for a given function / G L2{G) its interpolating 
spline Y^'^ is always an optimal approximation (modulo given information). 
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Theorem 2.4. The following statements hold true: 

1) IfK< {el + J'^ then the set Q{W, f, t, s, K) is empty. 

2) Every variational spline Y^'^ is the center of the convex set Q{W, f,t,e, K). 
As a result the following inequalities holds true for any g G Q{W, f, t, e, K) 



and 



w,f 



t,e 



9 



< -diamQ{W,f,t,e,K), 

HtAG) 2 



1 



(e/ + £)-*/2 diamQ{W,f,t,e,K), 



< 

L2(G) 2 

where diam is taken with respect to the norm of the Sobolev space Ht^e{G). 

Proof. Given a function / e L2{G) the linear manifold I{W,f) is the set of all 
functions g from L2{G) such that f{w) = g{w),w G W. Let us note that the 
distance from zero to the subspace I{W,f ), in the metric of the space Ht^e{G) is 
exactly the Sobolev norm of the unique spline Y^^J-^ e I{W, /). This norm can be 
expressed in terms of the sequence {Y^'-^){w) = f{w),w G W, and the sequence 
{aw{Y^'^)},w e W, from the representation 



yWJ 



Indeed, 



= ({ei+cr/^Y^/,iei+cr/^Yy'/ 



1/2 



{el + £)%':'f ,Y, 



1/2 




aUY,';'/)f{w) 



\wew J 
It shows that the intersection 

Q{W, f, t, e, K) = I{W, f) fl St,,(0, K), 

where -Bt,£(0, K) is the ball in Ht^s{G) whose center is zero and the radius is K, is 
not empty if and only if 



1/2 



K > 



Y 



WJ 



t.e 



HtAG) 



\wew J 



The first part of the Theorem is proved. 

Now wc arc going to show that for a given function / the interpolating spline 
Y^^^ is the center of the convex, closed and bounded set Q{W, f,t,e, K) for any 



K > \\Y 



t,e 



\Ht,e{G) ■ In other words it is sufficient to show that if 

Y,"^/ +heQ{WJ,t,e,K) 

for some function h from the Sobolev space Ht,e{G) then the function Y^'-^ — h 
also belongs to the same intersection. Indeed, since h is zero on the set W one has 
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(el + {si + Lf'h) = ({el + L)%'^/ , h) = 



But then 



i2(G) 



{el + Cy/^Y,"^/ + h) 
In other words, 

' {el + Lf^Y.'^'f - h) 



{el + Cy/^ {Y,^/-h 



L2{G) 



L-2{G) 



< K 



and because Y^'^ +h and Y^^'^ — h take the same values on W the function Y^'^ — h 
belongs to Q{W^ f, t, e, K). It is clear that the following inequality holds true 

\\yZ"^ - 9\\h,ag) < ldtamQ{WJ,t,e,K) 
for any g £ Q{W, /, t, e, K). Using this inequality one obtains 

{eI + LV\eI + Cf'' {YZ^f-g) 
(el + C)-*^^ diamQ{WJ,t,e,K). 



L2{G) 
1 



L2(G) 



< 



The Theorem is proven. 



□ 



3. Paley- Wiener spaces on combinatorial graphs 

The Paley- Wiener spaces PW^{G),lo > 0, were introduced in the Definition 1 
of the Introduction. Since the operator C is bounded it is clear that every function 
from L2{G) belongs to a certain Paley- Wiener space. Note that if 

Wmin = inf UJ 

then the space PWuj{G) is not trivial if and only if w > Wmin- 

Using the spectral resolution of identity P\ we define the unitary group of oper- 
ators by the formula 



The next theorem can be considered as a form of the Paley- Wiener theorem and 
it essentially follows from a more general result in |25) . 

Theorem 3.1. The following statements hold true: 

1) f £ PW^ {G) if and only if for all s £ the following Bernstein inequality 
takes place 

(3.1) < ^11/11; 

2) the norm of the operator C in the space PWui{G) is exactly uj; 

3) f £ PWi^{G) if and only if the following holds true 

lim \\£'f\\^/' =uj,s £R+; 
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4) f & PWuj (G) if and only if for every g <E L2 (G) the scalar-valued function of 
the real variable t 

{e'"'f,g) = Y.e'"'f{v)W) 

is bounded on the real line and has an extension to the complex plane as an entire 
function of the exponential type tu; 

5) f € PWui(G) if and only if the abstract-valued function e'*^/ is bounded on 
the real line and has an extension to the complex plane as an entire function of the 
exponential type uj; 

6) f € PWuj{G) if and only if the solution u{t,v),t G e V{G), of the 
Cauchy problem for the corresponding Schrodinger equation 

.du{t, v) 



- = Cu{t,v),u{0,v) = f{v),i = V-i, 

has analytic extension u{z,v) to the complex plane C as an entire function and 
satisfies the estimate 

ll"(^,-)IU.(G)<e"l^^l||/|U,(G). 
We prove here only the first part of the Theorem. 

Lemma 3.2. A function f e L2{G) belongs to PW„(G) if and only if the following 
Bernstein inequality holds true for all s e M+ 

(3.2) ||/:vil<c.^ll/l|. 



Proof. We use the spectral theorem for the operator C in the space L2{G) in the 
form it was presented in the Introduction. 

Let / belongs to the space PW^,(G) and J^cf = x G X. Then 

A2^||x(A)||5,(,)dm(A)) = (y^^ A2^||x(A)||5,(,)dm(A)) < u;^\\x\\x, s e M+, 

which gives Bernstein inequality for /. 

Conversely, if / satisfies Bernstein inequality then x = Tcf satisfies < 
a;**||a;||x- Suppose that there exists a set ct C [0, 00) \ [0, w] whose m-measure is not 
zero and x\a ^0. We can assume that cr C [w + e, 00) for some e > 0. Then for any 
s e K+ we have 

/ l|,x(A)|||(,)dm(A)< r A-2^||AMA)|l^(,)rfM<INIIWc^ + er, 
which shows that or a; (A) is zero on cr or a has measure zero. □ 

The Theorem 3.1 shows that the notion of Paley- Wiener functions of type u 
on a combinatorial graph can be completely understood in terms of familiar entire 
functions of exponential type w bounded on the real line. 

The notion of A-sets was introduced in the Definition 3 in the Introduction. The 
role of A-sets is explained in the following Theorem. 

Theorem 3.3. // a set S C V{G) is a A-set, then the set U = V{G)\S is a 
uniqueness set for any space PW^iG) with co < 1/A. 
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Proof. U f,g e PWu:{G) then f - g e PWu{G) and according to the Theorem 3.1 
the following Bernstein inequality holds true 

(3.3) \\C{f-g)\\<u;\\f-g\\. 

If / and g coincide on U = V{G)\S then f — g belongs to L2{S) and since 5 is a 
A-set we have 

||/-.g||<A||£(/-.9)||. 

Assume that w < 1/A and that / is not identical to g. We have the fohowing 
inequalities 

11/ - g\\ < A||£(/ - 5)11 < AcoWf - g\\ < \\f - g\\,co < 1/A, 

which provide the desired contradiction if / — is not identical zero. It proves the 
Theorem. □ 

As it was mentioned in the Introduction one cannot expect that non-trivial 
uniqueness sets there exist for functions from every Paley- Wiener subspace. But it 
is reasonable to expect that uniqueness sets exist for Paley- Wiener spaces PW^ (G) 
with relatively small u) > 0. Indeed, a direct calculation shows that for any graph 
G spaces PW^{G) with 



< w < Wl + — ^ = J^G > l,d(G) = max d{v), 
Y d{G) vev{G) 

have non-trivial uniqueness sets. 

Here are two examples of Paley- Wiener spaces on graphs and their uniqueness 
sets. 



1. Finite graphs. If a set of vertices V{G) of a graph G is finite then the 
spectrum of the Laplace operator is discrete and the space PWu{G) is a span of 

eigenfunctions whose eigenvalues < lo. In this case if ?7 is a uniqueness sets for a 
space PWu}{G) then \U\ is at least a number of eigenvalues (with multiplicities) of 
£ on the interval [0,uj]- 



2. Lattice Z". We consider a one-dimensional lattice Z. In this case there is 
a version of the Fourier transform on the space L'2{Z) which is defined by the 
formula 

Hfm = E f(^y'^' / e L2(Z), ^ e [-TT, tt). 
feez 

It gives a unitary operator from L2{G) on the space L2{T) = L2{T, where T 

is the one-dimensional torus and di^jl-K is the normalized measure. One can verify 
the following formula 

J-(/:z,/)(0 = 2sin2 \T{m). 

where £z is the Laplace operator on the graph Z. The next result is obvious. 

Theorem 3.4. The spectmm of the Laplace operator Cz on the one- dimensional 
lattice Z is the set [0,2]. A function f belongs to the space PWu(Z),0 < u) < 2, if 
and only if the support of Tf is a subset 0,^ of [— tt, tt) on which 2 sin^ 2—'^- 
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Our nearest goal is to show that for a one-dimensional line graph Z the esti- 
mates in Poincare inequalities of finite successive sets of vertices can be computed 
explicitly. 

Consider a set of successive vertices S = {vi,V2, vn} C Z, and the correspond- 
ing space L2{S). If bS = {vq, Vn+i} is the boundary of S, then for any ^p G L2{S) 
the function Cz'fi has support on S U bS and 

Czvivo) = -(p{vi),Cz(pivi) = 2ip{vi) - (p{v2), 
and for any other Vj with 2 < j < — 1, 

Cz(p{Vj) = -(p{vN-l) + 2(p{Vj) - (p(vN+l)- 

Let C2N+2 — ^{S) be a cycle graph 

C2N+2 = {U-N-1,U-N,---,U-1,Uo,Ui,U2,...,Un,Un+i} 

with the following identification 

U-N-l = Un+i. 

Thus the total number of vertices in C2N+2 is 2N + 2. We introduce an embedding 
of S U 65 into C2Ar+2 by the following identification 

ti(, = Uo,Vl = Ul, ...,Vn = Un,Vn+1 = Un+1- 

This embedding gives a rise to an embedding of L2{S) into L2(C2n+2), namely 
every € L2{S) is identified with a function € L2{C2n+2) for which 

F^{uo) = 0,F/^(wi) = ip{vi),...,F^{un) = fivN), F^{un+i) = 0, 

and also 

F^{u-i) = -ip{vi), ...,F^{u-n) = -ip{VN)- 
It is important to note that 

UGC2N+2 

If jCc is the Laplace operator on the cycle C2N+2 then a direct computation shows 
that for the vector defined above the following is true 

211^11 = ||F^||,2||£z(^|| = \\CcFj,^G L2{S),F^ e i2(C2iv+2). 

The operator jCc in L2{C2n+2) has a complete system of orthonormal eigenfunc- 

tions 

Tl 

(3.4) =exp27ri ^^_^^ A;,0< n < 2N + 1,1 < k < 2N + 2, 
with eigenvalues 

(3.5) A„ = l-cos^^^-^,0<n<2Ar + l. 

The definition of the function F^ € L2{C2n+2) implies that it is orthogonal to all 
constants and its Fourier series does not contain a term which corresponds to the 
index n = 0. It allows to obtain the following estimate 

2Ar+l 

\\CcF^r= E A^I(i^v',V'n)l'>4sin4^^||i^^f. 
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It gives the fohowing estimate for functions ip from L2{S) 

ll^ll<^sin-2FT2l|/:z^l|. 

Thus we proved the following Lemma. 

Lemma 3.5. If S = {vi,V2, ■■■,vn} consists of \S\ = N successive vertices of a 
line graph Z then it is a A-set for 

A = - sm 



2 2|S'| + 2' 

In other words, for any p> € L2{S) the following inequality holds true 

\M<M\^M- 

Remark 3. The last inequality which can be written as 

\\Czp\\ > 2 sin^ ^^^^M,p e L2{S), 

is similar to one of inequalities in [7] . 

Note that in the case |5| = 1 the last Lemma gives the inequality 

\M < ||«||,5 = M, 

but direct calculations give a better value for A: 

^\\C^S,\\,v&V. 

Let us note that if {Sj} is a finite or infinite sequence of disjoint subsets of 
vertices Sj C V such that the sets Sj U bSj are pairwise disjoint and every Sj has 
type Aj, then their union S = [J^ Sj is a set of type A = sup^ Aj. Indeed, since 
the sets Sj are disjoint every function p e L2{S), S = IJ^ Sj, is a sum of functions 
ipj G L2{Sj) which are pairwise orthogonal. Moreover because the sets Sj UbSj are 
disjoint the functions Cpj are also orthogonal. Thus we have 

3 3 

where A = sup^ Aj . 

A combination of this observation along with the last Lemma 3.6 gives the fol- 
lowing result for any < a; < ^3/2. 

Theorem 3.6. If S is a finite or infinite union of disjoint sets {Sj} of successive 
vertices such that 

1 ) the sets Sj ~ Sj U bSj are disjoint 
and 

2) for every j the following inequality holds 

(3-6) 1^,1 <- ^^-1, 

2 arcsm ^ ^ 

then every function f G PWtj(Z) is uniquely determined by its values on the set 
U = V{Z)\S. 
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A similar result holds true for a lattice Z" of any dimension. Consider for 
example the ease n = 2. In this situation the Fourier transform on the space 
1/2 (^^) is the unitary operator T which is defined by the formula 



^(/)(a,6)= T] /(fci,A;2)e*'=^«^+"=^«%/eL2(ZxZ), 



E 



where (^1,^2) G [— 7r,7r) x [— 7r,7r). The operator is isomorphism of the space 
L2{G) on the space ^2(1 x T) = L2{T x T, ^^1^^2/4772), where T is the one- 
dimensional torus, the following formula holds true 

H^z^fm = (sin' I + sin' |) J^ifM, 

where £^2 is the Laplace operator on the graph Z^. We have the following result. 

Theorem 3.7. The spectrum of the Laplace operator on the lattice I? is the set 
[0, 2] . A function f belongs to the space PW^{Z x Z), < w < 2, i/ and only if the 
support of J^f is a subset fl^ of [— 7r,7r) x [— 7r,7r) on which 

• 2 Ci , • 2 ^2 , 
sm — + sm — < CO. 
2 2 - 

Given a set S = {v^.m}, 1 < n < N,l < m < M, we consider embedding of S 
into two-dimensional discrete torus of the size T = {2N + 2) x (2M + 2) = 
Every / e L2{S) is identified with a function g G L2{T) in the following way 

g{Un,m) = f{Vn,m), 1 < 71 < N,l < m < M, 



and 

We have 



9{un,m) =0,N<n<N + 2,M<m<M + 2. 



where Ct is the combinatorial Laplacian on the discrete torus T. Since eigenfunc- 
tions of £t are products of the corresponding functions (3.1) a direct calculation 
gives the following inequality 

M<\ 7 ^||£z^(^IUeL2(5). 

min (^sin jj^, sin j 

In a similar way one can obtain corresponding results for a lattice Z" of any 
dimension. Note that the spectrum of the Laplace operator on Z" is [0,2] and 
Qzn = ^{2n + l)/2n. 

Let Nj = {Nij, ...,Nn.j},j € N, be a sequence n-tuplcs of natural numbers. For 
every j the notation S{Nj) will be used for a "rectangular solid" of "dimensions" 

Nij X N2J X ... X NnJ- 

Using these notations we formulate the following sampling Theorem. 

Theorem 3.8. // S is a finite or infinite union of rectangular solids {S{Nj)} of 
vertices of dimensions Nij x N2J x ... x N^j such that 

1) the sets Sj = S{Nj) U bS{Nj) are disjoint, 
and 

2) the following inequality holds true for all j 

A ■ ( ■ • 

w < 4min sm— — T,sin— — sm 



2Ari,,- + 2' 2Ar2j + 2' 2Ar„j + 2 
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then every f € PWij(Z") is uniquely determined by its values on U — V(1j"')\S. 
4. Reconstruction of Paley- Wiener spaces using splines 

Now we are going to use variational splines ^ as a reconstruction tool of 
Paley- Wiener functions / € PWuj{G) from their values on uniqueness sets of the 
form U ~ V{G) \ S, where S* is a A-set and A < 1/lu. We will need the following 
Lemma. 

Lemma 4.1. If A is a bounded self-adjoint positive definite operator in a Hilbert 
space H and for an if G H and a positive a > the following inequality holds true 

then for the same ip G H , and all k — 2\l — 0, 1, 2, ... the following inequality holds 

ll^ll <«'PVl|. 

Proof. By the spectral theory [3| there exist a direct integral of Hilbert spaces 

r\\A\\ 

X = / X{T)dm{T) 



and a unitary operator F from H onto X, which transforms domain of > 0, 
onto Xt = {x G X|r*a; G X} with norm 

II^Vlk= T''\\Ff{r)rx^^^dm{r)\ 
and = T*{Ff ). According to our assumption we have for a particular ip G H 

p\\A\\ 

\Fip{T)\^dm{T) <a'^ / T^\FLp{T)\^dm{T) 
Jo 

and then for the interval B = B{0,a^^) we have 

\F(f{T)\^dm{T) + I \F(f\^dm(T) < 

B -'[0,||A||]\S 

( [ T^\Fipfdm{T)+ [ T'^\Fipfdm{T)] . 

\Jb -'[o,||A||]\s J 

Since aV^ < 1 on B{Q,a-^) 

0< / {\Fipf - a'^r'^lFiff) dm{T) < [ {a^T^\Fip\'^ ~ \Fip\'^) dm{T). 

Jb J[oM\\]\b 

This inequality implies the inequality 

< / (aV2|F(^|2 _ aV|i^(^|2) dm{T) < f {a\^\F(p\^ - a^T^\F(p\^) dm^r) 

Jb J[oM\\]\b 

or 



'10,\\A\\] 

which means 

M < a\\Ap\\ < a'WAM- 
Now, by using induction one can finish the proof of the Lemma. The Lemma is 
proved. □ 
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Proof of the first part of the Theorem 1.1. 

We assume that the operator L has bounded inverse and e = 0. If / € PW^i (G) 
and Y;^'' , 

^k' —^ko}^^^ variational spline which interpolates / on a set {/ = 
V{G) \ S, where 5 is a A- set (0 < w < 1/A), then / - Y^'^ e 1.2(5) and we have 

(4.1) ||/-n^'^||<A||£(/-y,^'Ol|. 

U f 

At this point we can apply the last Lemma with A = C, a = A and = f — Yf.'. 
It gives the inequality 

(4.2) \\f-Y;^^f\\<K^\\L\f-Yl''^')\\ 

for all fc = 2', I = 0, 1, 2, ... Since the interpolant Y^'^ minimizes the norm ||£*^ • || 
it gives 

11/ - '■'^1! < 2A'=||/:^/l|, fc = 2Kie N. 

Because for functions / € PWi^{G) the Bernstein inequality holds 

||/:"/|| <c."||/||,meN, 
it implies the first part of the Theorem 1.1: 

11/ - if '^11 < 27i/||,7 = Aw < 1, fc = 2',Z e N. 
Proof of the second part of the Theorem 1.1. 

Now we assume that the operator £ is not invertible (it is a typical situation on 
any finite graph). We fix an 

< e < T 
A 

and assume that 

< w < — - £. 
A 

If / e PWu{G) and Yj/'J is a variational spline which interpolates / on a set 
U = V{G) \ S where 5 is a A- set then / - Y^^'/ e L^iS) and we have 

(4.3) ||/-y,^^||<A||£(/-lf/)||. 
For any g G L2{G) the following inequality holds true 

(4.4) \\Cg\\<\\{sI + C)g\\. 

Thus the inequalities (4.3) and (4.4) imply the inequality 
||/-n^^||<A||(e7 + £)(/-F^/)||. 



U f 

We apply the Lemma 4.1 with A = si + jC, a = A and 'P = f — Y^'J ■ 



gives 



the inequality 

\\f -Y;^J'\\< A^\\{el + C)\f - Yl;[f)\\ 

for all = 2', / = 0, 1, 2, ... Using the minimization property of Y^'^^ we obtain 

11/ - < 2A'=||(£/ + cffik = 2\ie n. 

li f G PWi^{G), then the Bernstein inequality 

||£"/||<a;"||/||,meN, 

implies the inequality 

\\{6l + Crf\\<{w + sr\\f\\,m€N. 



VARIATIONAL SPLINES AND PALEY- WIENER SPACES ON COMBINATORIAL GRAPHS 17 



After all we have the following inequality 

11/ - n"!/!! < 27'=||/||,7 - A(u; + £) < 1, fc = 2',/ € N. 
The proof of the Theorem 1.1 is complete. 
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